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ABSTRACT 


A Model Description Document for the Emulation Simulation Computer Model 
Dublished previously. The model consisted of a detailed model (emulation) of 
a SAWD CO? removal subsystem which operated with much less detailed 
( simulation^ models of a cabin, crew, and condensing and sensible hea 
exchangers. The purpose was to explore the utility of such an emulation/ 
simulation combination in the design, development, and test of a piece of AR 
hardware - SAWD. 

Fxtensions to this original effort are presented in the manual. The first 
extension 5 1 s an update' of the model to reflect changes In the SAWD control 
logic which resulted from test. In addition, slight changes were also made to 
the SAWD model to permit restarting and to improve the iteration technique. 
The second extension is the development of simulation models for more pieces 
of air and water processing equipment. Models are Presented for. E , 
Molecular Sieve, Bosch, Sabatier, a new condensing heat exchanger, SPE, SFWEb, 
Catalytic Oxidizer, and multifiltration. The third extension is to create two 
system simulations using these models. The first system presented consists of 

one air and one water processing system. The second system ’nJ!"! Jablt at 
potential Space Station air revitalization system complete with a habitat, 

laboratory, four modes, and two crews. 


' } 
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FOREWARD 


This Model Description Document has been prepared by Hamilton Standard 
Division of United Technologies Corporation for the National Aeronautics and 
Space Administration's Langley Research Center in accordance with Contract 
NAS1-17397 , "Development of an Emulation/Simulation Computer Model of a Space 
Station Environmental Control and Life Support System (ECLSS)". This manual 
describes the analytical models used in the three computer simulation programs 
developed under this contract. 

Appreciation is expressed to the Technical Monitors, Messrs. John B. Hall, Jr. 
and Lawrence F. Rowell of the NASA Langley Research Center for their guidance 
and advice. 

This manual was written by Dr. James L. Yanosy, Program Engineer, with 
assistance from Mr. Stephen A. Giangrande. The extensions to the program 
presented in this manual were performed under the direction of Mr. John M. 
Neel, Program Manager. Thanks is given to Mr. Joseph M. Homa for his efforts 
in the development of the Space Station Model. 
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1.0 INTRODUCTION 


The purpose of the original ESCM program was to demonstrate the 
utility of an emulation simulation computer program in the design, 
development and test of a piece of life support equipment. The piece 
of life support equipment selected for emulation was the SAWD CO 2 
removal subsystem. A continuation of the effort called for an update 
of the computer model following testing of the SAWD unit. In 
addition, extensions to the contract called for the development of 
"lightweight" or low fidelity simulation models of contending life 
support equipment and the configuring of this equipment into two 
different systems. 

This document provides three appendices to the original model 
description document[l]*. The appendices are: 

A. Emulation Simulation Computer Model Update 

B. Life Support System Model 

C. Space Station Model 

The following errata were found in the original document: 

(1) Cover page: Report number should be SVHSER 9504. 

(2) Page 25, fifth line: Should be "HA =* Total enthalpy of gas 

entering header, Btu" 


♦Numbers in brackets denote references listed in Section 2.0 


1 
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2.0 REFERENCES 


(1) Yanosy, J., "Model Description Document for a Computer Program 
for the Emulation/Simulation of a Space Station Environmental 
Control and Life Support System (ESCM)", Hamilton Standard Report 

SVHSER 9504 for NASA Langley Research Center, NASA CR- 181737, September 
1988. 

(2) "G189A Generalized Envi ronmental /Thermal Control and Life Support 
System Computer Program Manual", McDonnel Douglas Corporation 
MDAC-G2444; September, 1971. 

(3) Yanosy, 0., "User's Manual for a Computer Program for the 
Simulation of a Space Station Environmental Control and Life 
Support System (ECLSB) " , Hamilton Standard Report SVHSER 10630 
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A.l Introduction 

The Space Station Environmental Control and Life Support System 
(ECLSS) Emulation/Simulation Computer Model (ESCM) has been updated to 
include: (1) the capability of restarting the SAWD C0 2 removal 

subsystem from a transient start-up, (2) placing clamps on the SAWD 
bed segment temperature when iterating the bed segment temperature 
during the energy balance of the bed segment, and (3) the addition of 
an energy balance control method where the amount of C0 2 on the bed is 
determined from the time for the C0 2 to begin coming off the bed as 
detected by the flow sensor. This latter change was made to reflect 
a change made to the control logic in the hardware. Each of these 
updates to ESCM as they pertain to the Model Description Document are 

presented in this Appendix. 

A. 2 SAWD Bed 


The capability to start the SAWD subsystem without first running the 
steady state analysis was incorporated into the IR45 subroutine. The 
inlet and exit heater parameters, as well as SAWD bed segment 
parameters (i.e. temperature, pressure, molecular weight, and flow 
rates) were initialized in the IR45 subroutine; in a similar manner, 
they are initialized in the STEADY subroutine which is only called by 
IR45 in steady state. Therefore, the SAWD C0 2 removal subsystem can 
be started from a transient start-up as well as a steady state 

start-up. 


A-l 
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Additional logic was also added to the IR45 BALANCE subroutine. 
Limits were placed on the SAWD bed segment temperature while iterating 
on the bed segment temperature for an energy balance of the segment. 
One of the limits employed was as follows: if the temperature 

entering the bed segment was greater than the SAWD bed temperature at 

the start of the time step, then the resulting bed segment temperature 
must increase. 

A. 3 SAWD Control Model 

A new control method has been incorporated into the SAWD C0 2 removal 
subsystem model (IR45). This type of control uses an energy balance 
control scheme as opposed to the relative humidity method which as 
previously used. The energy balance method uses energy principles to 
determine the absorption time of the next cycle based on the bed's 
past desorb cycle. The total amount of energy for desorb can be 
calculated form the amount of steam added to the bed during desorb. 
The amount of C0 2 desorbed from the bed is known from techniques using 
the accumulator on the flow sensor; thus, the energy required to 
remove the C0 2 form the bed can be calculated. Also, the energy 
required to heat up the bed resins and canister before the desorption 
of C0 2 can be calculated. Therefore, from an energy balance, the 
amount of energy to heat the water which was on the bed at the start 
of the desorb can be calculated. Thus, the bed water loading at the 
start of the desorb can be determined. Knowing the amount of water on 
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the bed at the start of the desorb, the next absorption time can be 
determined. The following SAWD control logic has been incorporated in 
ESCM Subroutine GPOLY1 to determine the absorb time. 


Control Constants (Set in subroutine GP0LY1):] 


K0 - 1.11 
K1 - 0.045 
K16 - 115.01 
K41 - 0.25 


K42 - 2.5 
K43 - 60.0 
K44 - 3.808 
K46 - 0.6996 


K47 - 7480.0 
K48 - 60.87 
K49 - .00028205 
K50 * 0.0034246 


K51 - 1085 
K52 - 157.37 
K53 - 4.3167 
K54 - 0.040259 


Input: 

C02TIME 

DETIME 

INTI 

NEWABTIME 

ONTIME 

P 

PH20PAST 

TIN 

TOUT 


Time CO 2 goes to cabin during desorb, minutes 
Time for desorption, minutes 
Previous INTI value. 

Previous absorb cycle time, minutes. 

Time steam generator is on during desorption, sec. 
Accumulator pressure at end of desorption, sec. 

Past value of bed loading for this bed. %• 

Inlet temperature at end of desorbing bed's post absorb 

eye 1 e » F 

Exit temperature at end. of desorbing bed's past 
adsorbing cycle, F 
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All the following calculations are done at the end of desorption. 


Calculate Bed CO 2 loading fraction which was at start of this 
desorb "FC02". 


Power 

- 0NTIME/(0. 05769 * DETIME) 

Power 

=■ Clamp (Power, 1, 1040) 

TSAT 

- K54 * P**2 + K53 * P + K52 

T3 

- K 1 6 * (TSAT - 212) /Power 

T3 

- Clamp (T3, 0.0, 17.5) 

Ratio 

- (C02 Time - T3)/DETIME 

FC02 

- K51*RATIO**2+K50*RATIO-K49 

FC02 

- Clamp (FC02, 0.0, 0.05) 


Calculate Bed H 2 O loading percent which was at start of this desorb 
"PCTH20" . 

WGTH20 =■ (KO * ONTIME - K46 * (TSAT-TIN) - K47 * FC02 - K48)/ 
(TSAT-TOUT) - K1 * DETIME - K44 
0CTH20 - 100.0 * WGTH20 

IF( (PCTH20 - 0CTH20PAST) .LT.-3.0) PCTH20 * PCTH20PAST-3.0 
IF( (PCTH20-PCTH20PAST) .LT.-3.0) PCTH20 - PCTH20PAST-3.0 
PCTH20 - Clamp (PCTH20, 10.0, 40.0) 

Calculate new absorb cycle time, "NEWABTIME" in minutes. 
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INTI - K41 * (PCTH20 - 25.0) + INTI 

INTI - Clamp (INTI, -35.0, 35.0) 

NEWABTIME * K43 + K42 * (PCTH20 - 25.) 

NEWABTIME * NEWABTIME + INTI 
NEWABTIME - Clamp (NEWABTIME, 20., 90.) 


The steam flow rate determines the desorption time. The control of the 
steam flow is not calculated based on the relative humidity in the 
cabin during the previous absorption cycle. The steam flow for the 
next desorption of the bed is based on the steam flow rate during the 
past desorb multiplied by the calculated steam generation power ratio. 


Msn “ M SO * P R 

The steam generator power ratio is set equal to 1.0 for the first 
desorption of each SAWD bed. For all subsequent desorptions, the 
power ratio is calculated by squaring the actual time of the last 
desorption and dividing it by the product of the calculated time of 
the last desorption and the calculated time of the new desorption 

time. 


PR - fc a * fc a 


tco * ten 
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where: 


M$n = New bed desorb steam flow, pph 

M S0 = Old bed desorb steam flow, pph 

Pr * Steam generator power ratio 

fc a = Actual time of last bed desorb, sec 

t C o = Calculated time of last bed desorb, 

t cn ■ Calculated time of next bed desorb, 


sec 

sec 
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B.l Introduction 


An extension to the original ESCM program is to develop lightweight 
simulation models of various life support equipment and to combine 
them into a system. The utility of using all lightweight models to 
enhance system design could then be explored. This differs from 
the original program which was to investigate the utility of a 
combined emulation and simulation model. 


This manual provides the model description for one of the systems 
simulated. This system consists of one air revitalization group of 
equipment working in conjunction with one group of waste water 
processing equipment. The principal pieces of equipment are: 


Function 


Subsystem 


CO 2 Removal 
CO 2 Reduction 
O 2 Generation 
Trace Gas Removal 
Condensate Processing 
Urine Reclamation 


Electrochemical Depolarized Concentrator 
Sabatier 

Static Feed Solid Polymer Electrolysis 
Catalytic Oxidizer 
Multifi ltration 

Vapor Compression Distillation 


B-l 


SVHSER 10638 



UNITED 

TECHNOLOGIES 

H&MdlTW 


B.2 Model ling of System 


The life support system to be analyzed is shown in Figure B-l. The 
following discusses how this system represents real hardware. 


The system shown in Figure B-l consists of many components; 
however, some are the same. Therefore, models of only the 
following components are needed: 


(1) 

Crew 

(8) 

EDC Cells 

(2) 

Cabin 

(9) 

Fan 

(3) 

Fan 

(10) 

Splitter 

(4) 

Heat Exchanger 

(11) 

Mixer 

(5) 

SPE Cells 

(12) 

All Purpose Component 

(6) 

Catalytic Oxidizer 

(13) 

VCD 

(7) 

Sabatier 




These components are arranged as required for use with G189A. 
Accordingly, particular items are arranged for modelling purposes 
but do not represent their actual physical location. First of all, 
the crew which consumes oxygen and produces carbon dioxide and 
water vapor is placed in series with other components before the 
cabin. This was done to keep the schematic simple and to minimize 
the number of components required. In actuality, the crew is in 
the cabin. 
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In the model in Figure B-l, it appears that the cabin has only four 
ports for gas flows - two entry and two exit. This is a restric- 
tion of G189A and again does not represent actuality. Air for the 
condensing heat exchanger, Sabatier cooling and the EDC are all 
drawn directly from the cabin. Since G189A has a limited number of 
ports, one port is used and then splitters are used to direct flow to 
each of the components. The end result is the same. 


Modelling of Components 


As discussed in Section B.2, only thirteen component analytical 
models need be available to analyze the life support system. The 
following sections will discuss these analytical models, the 
generation of performance constants, and any other parameters 
necessary to describe the analytical model. In many instances, the 
analytical model is already described in the G189A manual [2]. In 
those cases, the reader will be referred to the G189A manual for a 
description of the analytical model. 

In addition to the thirteen components, the control of the cabin 


air conditions and the water tank levels are discussed 
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Of the thirteen components, many were discussed previously in the 
ESCM Model Description Document [1]. The following are discussed 

here: 

SPE Cells Bacteria Filter 

Catalytic Oxidizer Charcoal Filter 

Sabatier Multifiltration 

EDC Cells Controls 

VCD 

SPE Cells 


The solid polymer electrolysis cells convert water to hydrogen and 
oxygen and give off heat to a heat exchanger with the coolant being 
air. The following constants and operating conditions are used: 

A c “ unit cell area - 0.239 ft 2 
P c - cell operating pressure - 200 psia 
T c = cell operating temperature - 155°F 
N c * number of cells - 20 

The cell current density is determined from: 

Jq * I/A c 

where I - cell current, amps 

J c - cell current density, amps/ft 2 
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A cell efficiency is then determined by linear interpolation of 
the tables in Table B-l for the cell pressure P c , current density 
J c , and ceil temperature T c . in addition, the voltage V c across 
each cell is determined by linear interpolation of the tables in 
Table B-2 for the cell pressure P c , current density J c , and cell 
temperature T c . 


The total watts consumed by the cells becomes: 


w e - I V c N c 

where V c - voltage across each cell, volts 

w e - total electrical power of cells, watts 

The amount of oxygen and hydrogen produced is given by: 

m 0 2 - 0.000659 I N c 

T0<5 

m H2 * m 0 2/8 

where: mQ 2 - mass flow of oxygen produced, Ibm/hr 

m H2 * mass fl °w of hydrogen produced, lbm/hr 

The water consumed by electrolysis is given by: 

m H20,E M 1*125 mQ2 


B-6 
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TABLE B-l 

SPE CELL OPERATING EFFICIENCY 


(AMPS?FT 2 ) 

50 

100 

150 

200 

250 

300 

350 

400 

450 

500 


Percent Efficiency at T c - 

140°F 


Cell Pressure (psia) 



100 

200 

300 

91.1 

95.0 

96.8 

97.5 

98.0 
98.3 

98.6 

98.7 

98.9 

99.1 

84.7 

90.0 

93.5 

95.1 

96.0 

96.6 

97.2 

97.6 

97.8 

98.1 

82.7 

86.4 

90.3 

92.6 

94.0 
94.9 

95.7 
96.2 

96.7 

97.0 


Percent Efficiency at T c - 180°F 
Cell Pressure (psia) 


(AMPS?FT 2 ) 


50 

100 

150 

200 

250 

300 

350 

400 

450 

500 


100 


200 


300 


87.8 

93.0 
95.4 

96.6 

97.2 

97.7 

98.0 

98.3 

98.4 
98.0 


82.4 

86.7 

90.7 
92.9 
94.2 

95.1 

95.8 

96.4 
96.7 

97.1 


79.0 
82.9 
86.3 

89.1 

91.2 

92.6 

93.7 

94.5 
95.1 

95.6 
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TABLE B-2 
SPE CELL VOLTAGE 


Cell Voltage at T c - 140°F 
Cell Pressure (psia) 


(AMPS?FT 2 ) 

100 

200 

300 

50 

100 

150 

200 

250 

300 

350 

400 

1.556 

1.594 

1.623 

1.650 

1.678 

1.714 

1.760 

1.832 

1.568 

1.602 

1.632 

1.662 

1.690 

1.725 

1.772 

1.843 

1.577 

1.609 

1.638 

1.667 

1.697 

1.733 

1.777 

1.857 


Cel 1 Voltage at 

Tr - 180°F 
(psia) 



Cell Pressure 


(AMPS^FT 2 ) 

100 

200 

300 

50 

100 

150 

200 

250 

300 

350 

400 

1.503 

1.537 

1.562 

1.587 

1.614 

1.639 

1.670 

1.701 

1.514 

1.548 

1.573 

1.599 

1.624 

1.651 

1.682 

1.722 

1.522 

1.553 

1.581 

1.606 

1.632 

1.657 

1.689 

1.728 
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where: m^o, E “ mass flow of water consumed by electrolysis, 

lbm/hr 

Water is also present as saturated vapor in the oxygen and hydrogen 
gas flows. 

The amount is determined from the following relation: 


p sat( T c) 


' V - 

Pc 



m H20,02 “ 

Yv 

m 02 

MW h20 

1-Y V 

mw 02 


m H20,H2 * 

Yv 

m H2 

MW h2 

1-Y V 

MW H2 



where: Y v * mole fraction of water vapor 

MWq 2 » molecular weight of oxygen, lbm/mole 
MWh 2 ■ molecular weight of hydrogen, lbm/mole 
P sa t . saturation pressure at given temp., psia 
m H20,02 ■ mass flow °f water vapor In oxygen, lbm/hr 
m H20,H2 ■ mass flow °f water vapor In hydrogen, lbm/hr 

Therefore, the total water consumed becomes: 

m H20,T - m H20,E + m H20,02 + m H20,H2 
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The waste electrical power is given by: 
w w - w e (1 - 1.48 *i /V c ) 

The heat given off to the ambient is given by: 

Q - (w w + 0.11 w e + 27.3) 3.413 
where Q - heat lost to ambient, Btu/hr. 

All this heat is transferred to a heat exchanger which is cooled by 
ambient ai r. 

B.3.2 Catalytic Oxidizer 

The model of the catalytic oxidizer is a simple model which com- 
putes the temperature of the air stream exiting the unit. The 

composition of the air remains unchanged. The following constants 
are used: 

w htr * heater power - 28 watts 

M c « mass of high temperature bed catalyst - 2.0 lbm 
T c = operating temperature of high temp, catalyst * 600°F 
s - fraction of air flow to high temp, catalyst * 0.1111 
h L hx ■ heat exchanger effectiveness - 0.9 
Cp = specific heat of air - 0.24 Btu/lbm-F 
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First the temperature leaving the high temperature bed is computed 
as: 


Th - T c - ^ HX (Tc - T l) 

where - exit temp, from high temp, bed, °F 
T-j » air inlet temp, to cat ox., °F 

This high temperature air mixes with the air not going to the high 
temperature bed to give the resultant exit temperature. 

s mi T h c p + (1-S) mi Ti c p 

T e * 

mi Cp 

Figure B-2 shows a schematic of the catalytic oxidizer subsystem 
and the location of the above discussed temperatures and flows. 

B.3.3 Sabatier CO? Reduction Subsystem 

The Sabatier CO 2 reduction subsystem consists of the Sabatier 
reactor, a fan and a water separator as modeled for G189A. 

A listing of the Sabatier subroutine SABHS is given in the ECLSB 
User's Manual [3]. The model is a simplified block box model which 
provides the gas composition exiting the reactor and the mixed 
temperature of the cooling air leaving the reactor and condensing 
heat exchanger. 
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HTCO=High Temperature Catalytic Oxidizer 
ATCO=Ambient Temperature Catalytic Oxidizer 


FIGURE B-2 

CATALYTIC OXIDIZER SUBSYSTEM SCHEMATIC 
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The Sabatier reaction is assumed to go to completion. First a test 
is made on the incoming gases to determine whether excess hydrogen 
or carbon dioxide exists. The reaction is: 

CO 2 + 4 H 2 — *"CH4 + 2 H 2 O 
If m H2 ,i <4 m C 02,i 

2.016 44.01 

then excessive carbon dioxide exists and all the hydrogen is 
consumed. For this case, the following relations hold: 


m CH4,e * m H2,i MW CH4 
4 MW h2 

m H20,e * m H2,i 2 MW H20 + m H20,i 
4 MW H2 

m H2,e ■ °*° 

m C02,e * m C02,i “ m H2,1 MW C02 
4 MW H 2 

Qsen * 8864. * T *H2 , i 

4 MWh2 
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where: m^j = mass f1ow of hydrogen entering reactor, lbm/hr 

m H2,e ■ mass flow of hydrogen exiting reactor, lbm/hr 
m C02,i ■ mass flow of carbon dioxide entering reactor, 
lbm/hr 

m C02 , e ■ mass flow of carbon dioxide exiting reactor, 
lbm/hr 

m CH4,e ■ mass flow of methane exiting reactor, lbm/hr 
m H20,e ■ mass flow of water exiting reactor, lbm/hr 
m H20,i * mass flow of water entering reactor, lbm/hr 
MW » molecular weight, lbm/mole 
Qsens ■ sensible heat produced by reaction, Btu/hr 


On the other hand, if excess hydrogen exists, the following 
relations hold: 


m CH4,e - m C0 2,i MW CH4 

MWco2 

m H20e * m C0 2,i 2MW H2 o + m H2 o,i 
MW c0 2 

m C02,e * 0-0 

m H2,e “ m H2 ,i - mco2,1 4MW H2 
MW c02 
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Qsens 3 8865 * m C02 , i 4MW H2 

MWco2 

The latent heat load to condense the water produced is 

Q lat - 1050. m H2 o,e 

The temperature of the cooling air leaving the Sabatier is: 

T air,e - T a1r,i + (Qlat + Qsens) /m air c p 
where: m a i r - mass flow of air, lbm/hr 

The fan draws air around the reactor for cooling and through a 
condensing heat exchanger to condense the product water. The fan 
operational characteristics are: 

cfm - 20.0 
power » 34 watts 

The condensed water and product gases pass through a water 
separator which is modeled as an alternate component for use with 
G189A. All the water is separated from the product gases and 0.13 
Btu/hr is added to the product gases in the form of heat. 
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B.3.4 EDC CO? Removal Subsystem 


The EDC C0 2 removal subsystem removes carbon dioxide from the air 
stream through a fuel cell type process which requires hydrogen and 
produces electricity. Heat Is removed through a water cooled heat 
exchanger. Again, simple black box models are used. See Figure 
B-3 for a schematic representation. 


A listing of the EDC subroutine EDCHS is contained in Appendix A of 
the ECLSB User's Manual [3]. The following constants are used: 

0d * design current density » ll Amps/ft 2 

r C 02 m co 2 transfer rate - 0.001736 lbm C0 2 /amp-hr 

A c * area per cell » ft 2 /cell 

T c - cell operating temperature - 70°F 

P c = cell operating pressure * 14.7 psia 

c PC02 * 0*211 Btu/lbm-F 

C PH2 * 3.437 Btu/lbm-F 

CpQ2 = 0.221 Btu/lbm-F 

Nc - number of cells * 30 

First, the carbon dioxide removal at design conditions is 
calculated using the following: 

m C02,d - r C 02 J d A c N c 
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Then, a factor is computed for use in determining the actual CO 2 
removed. This factor is given by: 


J C02 



The actual CO 2 removed becomes: 


m C02,A * m C02,d k C 02 
m C02 , e “ m C02 , i " m C02,A 

where: J A - actual current density, amps/ft^ 

m C02,i “ mass flow °f CO 2 EDC, lbm/hr 
m C02,e ■ mass flow of CO 2 exiting EDC in air stream, lbm/hr 


The actual electrical characteristics required are: 


I * 1.27 Kco2 J D A C 
V - J A 0.5 N c 

dividing line on separate half space 


J D 

w - I V 
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In the process, oxygen is consumed; the amount is given by the 
following relations: 


PC02 " 


m C02,i p i 760 
MWco2 14.696 


7 


1.02 + 0.19 

PC02 


0.718 


2 

P C02 


MC02,C 


m C02,A MW 02 

2 MW C0 2 


m 02,e m m 02,i " m 02,C 


where: 


p. > inlet pressure of air to EDC, psia 

i"C02,i ” in1et mass flow of C ° 2 int0 EDC ’ lbm/hr 
MWco2 * molecular weight of CO 2 ■ 44.01 lbm/mole 

P C 02 = partial pressure of CO 2 in inlet air, mm Hg 
- effectiveness of conversion 
MWq 2 * molecular weight of oxygen * 32 lbm/mole 
m 0 2 c - mass flow of oxygen consumed, lbm/hr 
mo2 ,i * mass flow of oxygen into EDC, lbm/hr 
("02 e - mass flow of exiting EDC air stream, lbm/hr 
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The factor of 2 in the preceding equation arises from the overall 
reaction taking place in the EDC cells: 


CO 2 + 1/2 O 2 + H 2 > CO 2 + H 2 O + electricity + heat 


Accordingly, the amount of hydrogen consumed is given by: 


m 02c m Wh2 

m H2,C “ 2 

MW 02 

The amount of water produced is: 

m m H2,c mw H20 

m H20,P - 

MW H2 

Some of the water available leaves with the hydrogen and carbon 
dioxide gases. It Is assumed that the vapor is saturated and that 
the pressure in the H 2 -C0 2 exit port is 20 psia. The amount of 
water leaving with the H 2 -C0 2 mixture is given by the following: 

m H2,e * m H2 ,i ~ m H2,C 

mWmix - mH2,e + m C02,A 
m H2,e + m C02 
mw H2 MWqq2 
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Py * p sat( T c) 


MW H 20 

MWmix 


P V 

P T - P V 


Mh 20,VE = OJ l m C02,a + m H2,e) 

m H2 0,e - m H20,p + m H20,i + m H20,Vi " m H20,ve 


where: 


m H2,i 

m H2,e 

MW H2 

MW H2 0 

MWmix 


p sat 

cjj 

PT 

m H20,ve 

m H20,i 

m H20,vi 


mass flow of hydrogen into EDC, Ibm/hr 
mass flow of hydrogen exiting EDC, lbm/hr 
molecular weight of hydrogen, lbm/moles 
molecular weight of water, lbm/moles 
, molecular weight of H2-CO2 mix exiting EDC, 
lbm/moles 

■ vapor pressure of water in H2-CO2 stream, psi 

* saturation pressure for given temp, psia 
» humidity ratio 

* total pressure in H2-CO2 exit line, psia 

■ vapor leaving with H2 end C02» lbm/hr 

- mass flow of water into EDC from air stream, 

lbm/hr 

* mass flow of water into EDC with H2 stream, 


lbm/hr 
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The total heat generated is given by: 


Qd “ 65.7 + 766.8 niC02,d 

Qsens * (1.06 <C02 ~ 0.06) Q d 

where: Qd - design heat production, Btu/hr 

Qsens m actual heat production, Btu/hr 

Of this heat, 25% is added to the air stream and 75% must be 
removed by the heat exchanger. Accordingly, the exit air 

temperature is given by: 


0.25 Qsens 

T e ■ T 1 

( m air,i + m H 20,vi) ep 

where: » inlet air temperature, F 

m air,i * inlet air flow, lbm/hr 
c p =■ specific heat of inlet air - 0.24 Btu/lbm-F 

VCD Water Processing Subsystem 

Vapor compression distillation process purifies water through a 
distillation type process. The model at present is simple and 
determines the amount of power consumed and the amounts of water 

and brine delivered. The following relations and constants are 
used: 
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m H20,e “ °- 96 m H20 , i 


mb - 0.04 mn20,1 


where: 


w 


m H20,e 

40.1 


120 


m H20 i * mass flow of water into VCD - 2.30 lbm/hr for 3 
man units or 2.73 lbm/hr for 6 man unit. 

mH2 o, e * mass flow of c1ean water exit1ng VCD ’ lbm/hr 
m b - mass flow of brine, lbm/hr 
w - power required, watts 


B.3.6 Filtration Models 

Three filtration devices are modeled in ECLSB by the use of the 
Alternate Component subroutine. The three components are the 
bacteria filter, charcoal filter, and multifiltration. For each of 
these, the inlet flows and conditions pass through unchanged. 


B.3.7 Control 


A variety of items require control to operate the system 
effectively. C0 2 and 0 2 levels in the cabin atmosphere need to be 
maintained, and water levels in the various tanks need to be 
regulated. This control logic is contained in subroutines GP0LY1 
and GP0LY2, and listings of these are contained in Table B-3 and 
Table B-4 respectively. 
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TABLE B-3 

LISTING OF GPOLY1 FOR ECLSB MODEL 

SUBROUTINE GPOLY1 
C 
C 

C THIS SUBROUTINE PROVIDES THE CONTROL LOGIC FOR SYSTEM «B” ECLS 

C SIMULATION USING THE G189A COMPUTER PROGRAM. 

C 

C 

C TYPE STATEMENTS: 

C 

INTEGER ONCE 
REAL MR.INOM 

LOGICAL STEADY , CYCLIC , LTSIDE , OPEN , VCDON , MFTON 
DIMENSION STATEMENTS: 

DIMENSION V( 1 ) ,K( 1 ) 

DIMENSION TURINB ( 8 > , THWASH C 1 1 > , TSHOWR ( 2 ) . TDRINK ( 8 ) , TFOODP ( 5 ) 
COMMON STATEMENTS: 

COMMON /COMP/ DS( 15) ,N,NA1 ,NB1,NC,NCAB,NCFL,NEXT,NEXV,NK, 

1 NKEX , NKS , NKT , NLFL , NP , NPASS , NPF , NPFK6 ) , NO . NS , NSF , NSFT ( 6 ) , 

2 NSTR( 18) , NSUBR.NV.NVT, Y( 12) 

COMMON /RARRAY/ IMAXR,R(0250) 

COMMON /ECLST1/ KCH0UT,KPRNT.KPTINV(4>,KWIT,KWIT1,KWIT2, 

1 KWIT8 , KWIT4 , NUFF , KSTEDY 
COMMON /KANDV/ R 

COMMON /MISC/ DTIME,GRAV,KFLSYS,KOUTPT,KPDROP,KSYPAS,KTRANS, 

1 LPSUM ( 5 > , MAXCI , MAXLP , MAXSLP , MAXSSI , NCOMPS , NEWDT , NLAST , NPASPD 

2 MINSSI.PCMIN.PLMIN. START, STEADY. TIME. TIMEMX.TMAX.TMIN.WTMAX 
COMMON /CASE/ NCASE.NRSCS.NRECS 

COMMON /PROPTY/ CP0,CP(99> ,CPC0NL t CPC0NV,CPC02.CPDIL f CPOXY.CPTC 

1 GAMGAS , RHOO , RHO ( 99 > , VISCO , VISC ( 99 > , VISGAS , WTMO , WTM ( 99 ) , WTMCON 

2 WTMDIL , WTMTC , XKO , XK ( 99 ) , XKGAS , XKLIQ , VISLIQ 
COMMON /SOURCE/ A( 19) , B( 19 > ,CPA ,CPB, IA1, IB1 ,NA .NB.NPFS, NPFSTfG > 

1 NSFS , NSFSTf 6 ) , RHOA , RHOB , VISCA , VISCB , WTMA , WTMB , XKA , XKB 
COMMON /VLOC/ IP, IS, IC , IQ, IV , IVT, IEX, INEXK 
COMMON /LRC/ IDATE(2 ) , ISCHM 
C 

C DATA INITIALIZATION: 

C 

DATA TURINE / 0.1,8.0,6.0,9.0,12.0,15.0.20.0,25.0/ 

DATA THWASH / 0. 1 ,8. 0,5. 0,6.0, 9.0,10.0,11.0,12.0,15.0,20.0,25 / 
DATA TSHOWR / 15.0,25./ 

DATA TDRINK / 0.1,8.0,6.0,9.0,12.0,15.0,20.0,25.0/ 

DATA TFOODP / 0.1, 5., 10., 14., 26./ 

DATA KU.KH.KS.KD.KF / 1, 1, 1, 1, 1/ 

C 

EQUIVALENCE < V ( 1 > , K ( 1 ) > 
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INITIALIZE AT START OF STEAD! STATE SOLUTION. DONE ONLY ONCE. 

IF (STEADY .AND. KSYPAS .LE. 1) THEN 
LTSIDE - .TRUE. 

TCABO - 0.0 
TDESO - 0.0 
IFREQ - VV(2, 184) 

NMEN - KK( 1 , 16) 

H20ADD - VV(2, 128) 

CYCLIC - VV(2, 185) .GT. 0.9 
END IF 

IF(N .EQ. 1) THEN 

READ TABLE OF METABOLIC RATE VS MISSION TIME (24 HR CYCLE). 


TIMCYC 

- 

AMOD (TIME, 86400. > 

MR 

- 

VALUE (1, TIMCYC, 0.0) 

TCAB 

- 

VV(2, 104) 

TCAB 

- 

70. 

Q L 

- 

MR-480. ♦( MR/ 1000. + 10. )*<TCAB-60.0> 

QLMIN 

- 

0 . 22*MR+2 . 6* ( TCAB-60 . 0 ) 

QL 

- 

AMAXl(QL.QLMIN) 

QS 


MR-QL 

R( 66 ) 

- 

QS 

R(67) 

- 

QL 

R(82> 
END IF 


MR 

IF(N .EQ. 

2) 

THEN 


IF (STEADY .OR. MOD( KSYPAS. IFREQ) .EQ. 0) ISCHM - 1 
R ( 181) - DTIME 
R( 182) - DTIME/60. 


LIGHTSIDE - DARKSIDE? 

IF (CYCLIC) THEN 

TIMORB - AMOD (TIME, 6400. ) 
LTSIDE - TIMORB .LE. 2700. 
END IF 


NITROGEN ADDITION RATE 

R( 166) - 0.0 

WCN - 0.0 

PT - R(4) 

P02 - R(94) 

IF(PT .GE. 14.819) GO TO 220 

IF (PT .GE. 14.818 .AND. P02 .GE. 8.28) GO TO 220 
IF ( P02 .LT. 3.09) GO TO 220 
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REGULATOR LOGIC. 

IF (OPEN) GO TO 210 

N2 OPENING CURVES. 

WCN - VALUE ( 22 , PT , 0 . 0 ) 

OPEN - .TRUE. 

GO TO 220 

N2 CLOSING CURVES. 

WCN - VALUE ( 28 , PT , 0 . 0 > 

IF (WCN .LE. 0.0) OPEN - .FALSE. 

R( 166) - WCN 

OXYGEN AND H20 VAPOR ADDITION FROM SPE, LBM/HR 

R( 128) - H20ADD+VV (20,68) 

R( 166) - V( IV*2) 

TRACE CONTAMINANTS ADDITION 


C 


C 


c 


END IF 

IF(N . EQ. 8) THEN 

IF (.NOT. STEADY .OR. KSYPAS .GE. 4) THEN 


WSEC 

- VV ( 7 , 1 ) 

WPRI 

- VV(5, 1) 

WTOT 

- WSEC+WPRI 

SR 

- WSEC/WTOT 

R( 65 ) 

- SR 

IF 



END IF 

IF(N .EQ. 6) THEN 
R(84) - 0.0 

IF( VV(7,G8) .GE. B60.) R(84) - 1.0 
END IF 

IF(N .EQ. 9) THEN 

IF (.NOT. STEADY .OR. KSYPAS .GE. 4) THEN 
WSEC - VV( 14 , 1 ) 

WPRI - VV( 16, 1 )*VV( 28 , 1 )+VV(21 , 1 ) 

WTOT - WSEC+WPRI 

SR - WSEC/WTOT 

R(65) - SR 

END IF 
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END IF 
C 

IF(N .EQ. 11) THEN 

W8 - A( 1 )*R(72)/CPA 

W8 - AMAXKW8, A( 1 ) ) 

R(66> - EXP(2. 809-0. 00169*B( 1) )*W8*(0. 52+0. 00026*B( 1) ) 

END IF 

EXTRACT CONDENSATE WATER. 

IF(N .EQ. 18) THEN 
R(67) - A(7 ) 

A( 1 ) - A(l)-A(7> 

A ( 7 > - 0.0 

CPA - (A(5)*A(8)+A(6)*CPC0NV)/A( 1) 

END IF 
C 

IF(N .EQ. 15) THEN 

IF (.NOT. STEADY .OR. KSYPAS .GE. 4) THEN 
WSEC - VV( 16, 1) 

WPRI - VV(28, 1 )+VV(21 , 1 ) 

WTOT - WSEC+WPRI 

SR - WSEC/WTOT 

R(65) - SR 

END IF 
END IF 
C 

IF(N .EQ. 16) THEN 
R<84) - 1.0 

IF( .NOT. LTSIDE) R(84) - 0.0 
END IF 
C 

17 IF(N .NE. 17) GO TO 1799 

IF( STEADY .AND. KSYPAS .LT. 4) GO TO 1799 
WSEC - VV(21 , 1 ) 

WPRI - VV(28 , 1 ) 

WTOT - WSEC+WPRI 
SR - WSEC/WTOT 
R(65) - SR 

1799 CONTINUE 
C 

18 IF(N .NE. 18) GO TO 1899 
R(84) - 1.0 

IF( .NOT. LTSIDE) R(84) - 0.0 
1899 CONTINUE 
C 

19 IF(N .NE. 19) GO TO 1999 
R(84) - 1.0 

IF( .NOT. LTSIDE) R(84> - 0.0 
1999 CONTINUE 
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20 IF(N .NE. 20) GO TO 2099 
INOM - R ( 72 > 

P02 - VV(2,94) 

R(69 > - INOM 

IF ( P02 .GE. 8.1) R(69) - 0.9*IN0M 
IF ( P02 .LE. 2.9) R(69) - 1.1*IN0M 
CALL SK(1,20, 16) 

IF (LIS IDE) GO TO 2099 
R(69 ) - 0.0 
CALL SRC 0,20, 16) 

2099 CONTINUE 
C 

21 IF(N .NE. 21) GO TO 2199 
R(68> - VV(2Q,6S) 

2199 CONTINUE 
C 

28 IF(N .NE. 28) GO TO 2899 
PC02 - VV(2, 100) 

R(68> - R( 69) 

IF ( PC02 .GT. 8.0) R(68) - R ( 69 > 
IFCPC02 .LT. 2.0) R(68) - R(69) 
IF(LTSIDE) GO TO 2899 
R(68) - 0.0 
2899 CONTINUE 
C 

24 IF(N .NE. 24) GO TO 2499 
R(6S) - VV(28,65) 

2499 CONTINUE 
C 

28 IF(N .NE. 28) GO TO 2899 

EXTRACT CONDENSATE WATER. 

R(67) - A(7 ) 

A(l) - A( 1 )-A(7) 

A(7> - 0.0 

CPA - ( A ( S ) *A ( 8 > +A ( 6 ) *CPCONV >/A ( 1 ) 
2899 CONTINUE 
C 

82 IF(N .NE. 82) GO TO 8299 
CPA - CPCONL 
RHOA - RHO(l) 

CPB - CPCONL 
RHOB - RHO(l) 

WTMA - WTM(l) 

WTMB - WTM(l) 

VISCA - VISCd) 

VISCB - VISCd) 

XKA - XK( 1 ) 
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XKB - XK( 1 ) 

8299 CONTINUE 
C 

85 IF(N ,NE. 85) GO TO 8599 

URINE DUMP TO STORAGE TANK, LBM/HR 

A ( 1 ) - 0.0 

IF (NMEN .EQ. 8) WDMP - 17. /7. 

IF (NMEN .EQ. 6) WDMP - 26.28/7. 

TIMCYC - AMOD (TIME, 86400. > 

IFtTIMCYC .LT. TURINE(7)*8600. .AND. KU .EQ. 8) KU - 1 
IF (TIMCYC .LT. TURINE(KU>*8600. > GO TO 8510 
KU - KU+1 

KU - MIN0(KU,8> 

A( 1 ) - WDMP/DTIME*8600 . 

WASTE HAND WASH WATER DUMP TO STORAGE TANK. LBM/HR 
WHWASH - 0.0 

IF (TIMCYC .LT. THWASHt 10 >*8600. .AND. KH .EQ. 11) KH - 1 
IFtTIMCYC .LT. THWASH(KH)*8600. ) GO TO 8620 
KH - KH* 1 

KH - MINO(KH.ll) 

WHWASH - 11.5/10. 

A( 1 ) - A ( 1 > ♦WRWASH/DTIME*8600 . 

WASTE SHOWER WATER DUMP TO STORAGE TANK, LBM/HR 

8520 WSHOWR - 0.0 

IFtTIMCYC .LT. TSHOWRt 1 )*8600. .AND. KS .EQ. 2) KS - 1 
IFtTIMCYC .LT. TSH0WR(KS)*8600. ) GO TO 8580 
KS - KS*1 

KS - MIN0(KS,2> 

WSHOWR - 22.6 

A(l> - A(1)+WSH0WR/DTIME*8600. 

FLOW EXITING URINE WASH WATER STORAGE TANK, LBM/HR 

8580 R( 1 ) - 0.0 

IF(R(69) .LE. 80.) GO TO 8540 
VCDON - .TRUE. 

GO TO 8550 

8540 IF(R(69) .GE. 27.) GO TO 8550 

VCDON - .FALSE. 

8550 IF (.NOT. VCDON) GO TO 8599 
IF (NMEN .EQ. 8) R(l) - 2.80 
IF (NMEN .EQ. 6) R(l) - 2.78 
8599 CONTINUE 
C 
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41 


4199 

C 

42 


4299 

C 

46 


4610 


4699 

C 

66 


6699 

C 

68 


IF(N .NE. 41) GO TO 4199 
H20SPE - VV(20,67) 

R(l> - ( WHWASH+WSHOWR )/DTIME*8600 . +H20SPE+WMFLT 
CONTINUE 

IF(N .NE. 42) GO TO 4299 

WSEC - ( H20SPE+WMFLT ) *DTIME/8600 . 

WTOT - WSEC ♦WHWASH+WSHOWR 
IF ( WTOT .LE. 0.0) GO TO 4299 
SR - WSEC/WTOT 
R<65> - SR 
CONTINUE 

IF(N .NE. 46) GO TO 4699 
R( 1 > - 0.0 

WDRINK - 0.0 

TIMCYC - AMOD (TIME, 86400. ) 

IF ( TIMCYC .LI. TDRINK(7)*8600. .AND. ED .EQ. 8) KD - 1 
IF (TIMCYC .LT. TDRINK(KD)*8600. > GO TO 4610 
IF (NMEN .EQ. 8) WDRINK - 19.56/7. 

IF (NMEN .EQ. 6) WDRINK - 81.82/7. 

KD - KD+1 

KD - MIN0(KD,8> 

R( 1 ) - WDRINK/DTIME*8600 . 

FOOD PREPARATION WATER DUMP TO STORAGE TANK, LBM/HR 
WFOODP - 0.0 

IF (TIMCYC .LT. TFOODP ( 4 ) *8600 . .AND. KF .EQ. 6) KF - 1 
IF (TIMCYC .LT. TFOODP (KF) *8600. ) GO TO 4699 
KF - KF+1 

KF - MIN0(KF,8> 

IF (NMEN .EQ. 8) WFOODP - 4.74/4. 

IF (NMEN .EQ. 6) WFOODP - 9.48/4. 

R( 1 ) - R( 1 )+WF00DP/DTIME*3600. 

CONTINUE 

IF(N .NE. 56) GO TO 5699 

WSEC - WMFLT 

WTOT - WSEC+H20SPE 

IF ( WTOT .LE. 0.0) GO TO 5699 

SR - WSEC/WTOT 

R(65) - SR 

CONTINUE 

IF(N .NE. 58) GO TO 5899 
WMFLT - 0.0 
R( 1 > - 0.0 

WPOT - VV(46,69) 

IF(WPOT .GE. 40.) GO TO 5810 
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MFTON - .TRUE. 

CO TO 5S20 

5810 IF ( WPOT .LE. 45.) GO TO 5820 
MFTON - .FALSE. 

5820 IF (.NOT. MFTON) GO TO 6899 

IF ( R ( 100 ) *DTIME/8600 . .LE. R(69)> GO TO 6880 
R(l) - 0.0 
1MFLT - R(100) 

GO TO 5899 
5880 R ( 1 ) - R(100) 

5899 CONTINUE 
C 

61 IF(N .NE. 61) GO TO 6199 

IF (STEADY .AND. KSYPAS .LT. 4) GO TO 6199 
TCAB - VV(2, 104) 

IDES - V( IV+28) 

TTOL - 0.1 

IF( ABS ( TCAB-TDES ) .LT. TTOL) GO TO 6199 
A1 - 0.025 
ITER - 1 

C ALL ESTIM ( R ( 66 ) . TC AB , TDES , R650LD , TC ABO , TDESO . A 1 . ITER , NSTR ( 1 ) ) 
R(65) - AMAX1(R(66) ,0.0) 

R(66> - AMIN1(R(65) ,0.9) 

6199 CONTINUE 
C 

RETURN 

END 
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TABLE B-4 


C 


c 

c 

c 

c 

1 


LISTING OF GP0LY2 FOR ECLSB MODEL 
SUBROUTINE GP0LY2 

COMMON /COMP/ DS ( 15 ) , N , NA 1 , NB 1 , NC , NCAB , NCFL , NEXT , NEXV t NK , 

1 NKEX , NKS.NKT, NLFL , NP , NPASS , NPF . NPFT ( 6 ) , NQ , NS , NSF , NSFTt 6 ) . 

2 NSTR( 18) ,NSUBR,NV P NVT, Y( 12) 

COMMON /RARRAY/ IMAXR,R(0250) 

COMMON /ECLST1/ KCHOUT , KPRNT , KPTINV ( 4 > , KWIT , KWII 1 , KWIT2 , 

1 KWIT3 , KWIT4 , NUFF , KSTEDY 
COMMON /KANDV/ K 

COMMON /MISC/ DTIME , GRAV , KFLSYS , KOUTPT , KPDROP , KSYPAS , KTRANS , 

1 LPSUM ( 5 ) , MAXCI . MAXLP , MAXSLP , MAXSSI , NCOMPS , NEWDT , NLAST , NPASPD , 

2 MINSSI , PGMIN , PLMIN , START , STEADY , TIME , TIMEMX , TMAX , TMIN, WTMAX 
COMMON /CASE/ NCASE , NRSCS , NRECS 

COMMON /PROPTY/ CP0,CP(99> t CPC0NL,CPC0NV,CPC02,CPDIL,CP0XY.CPTC, 

1 GAMGAS , RHOO , RHO ( 99 > , VISCO , VISC ( 99 > , YISGAS . WTMO , WTM ( 99 ) , WTMCON , 

2 WTMDIL , WTMTC , XKO , XK ( 99 > , XKGAS , XKLIQ , VISL IQ 

COMMON /SOURCE/ A( 19) ,B( 19 > .CPA.CPB, IA1 , IB1 , NA.NB ,NPFS, NPFSTt 6) , 
1 NSFS , NSFST ( 6 ) , RHOA , RHOB , VISCA . VISCB , WTMA . WTMB , XKA , XKB 
COMMON /VLOC/ IP, IS, IC , IQ, IV, IVT, IEX.INEXK 

DIMENSION V( 1 ) ,K( 1 ) 

EQUIVALENCE (V(1),K(1>) 

LOGICAL STEADY 

IF(N.NE.l) GO TO 199 


CALC NET FLOWS DUE TO CABIN PRI & SEC FLOW LOOPS 

SBC02 - 0.0 
SBH20 - 0.0 


CALC NET H20 VAPOR CHANGE 


125 XH20 - R(70) - SBH20 - VV( 18,67) ♦ VV(28,75) ♦ VV(20,68> 

CALL SV(XH20,2, 187) 

CALC NET 02 CHANGE 

X02 - -R(68) - VV(28,78) VV(20,66> 

CALL SV(X02,2, 175) 

CALC NET C02 CHANGE 

XC02 - R(69 ) - SBC02 - VV(28,79) 

CALL SV(XC02,2, 177) 

1 99 CONTINUE 

C 

2 IF (N .NE. 2) GO TO 299 


B-25 



SVHSER 10638 


R(2) - 70. 

R(21) - 70. 

R(89> - 0.42 

R<98> - 46. 

R( 104 > - 70. 

299 CONTINUE 
C 

18 IF(N.NE. 18) GO TO 1899 
R(20> - R<67> 

R(21> - R 1 2 > 

R(22) - 24.7 
R(23) - 24.7 
1899 CONTINUE 
C 

28 IFtN.NE. 28) GO TO 2899 
R(20) - R(67) 

R ( 2 1 ) - R<2> 

R(22) - 24.7 
R(28) - 24.7 
R(68) - VV(26,68) 

2899 CONTINUE 
C 

RETURN 

END 
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B.3.7.1 Nitrogen Addition Control 

Nitrogen is added to maintain the total pressure in the cabin at 
a level of 14.813 psia. The same control as used and described in 
the ESCM Model Description Document [1] is used in the ECLSB model. 

B.3.7.2 Oxygen Production Control 



The oxygen level is controlled by adjusting the electrical current 
to the SPE cells according to the following: 

p 02 > 3.1 Psia I - 0.9 I N0M 

p 02 < 2.9 Psia I - 1.1 l N0M 

where: Pq 2 =■ partial pressure of O2 in cabin, psia 

I = current to SPE cells, amps 
r NOM “ nominal current to cells - 22 amps 

This nominal electrical current corresponds to a nominal oxygen 
consumption rate for three men at 0.255 lbm/hr. 
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B.3.7.3 CO? Partial Pressure Control 


At present, the EDC production current is not adjusted in response 
to changes in the CO 2 partial pressure in the cabin. The current 
density is maintained at a constant 11.0 amps per square foot. The 
user may change this control logic by making desired changes in the 
GP0LY1 subroutine. 


B.3.7.4 Cabin Temperature and Humidity Control 

Control of cabin temperature and humidity is accomplished by 
regulating the fraction of air flow from the cabin which passes 
through the condensing heat exchanger. The more air through the 
heat exchanger, the cooler the cabin should become. The technique 
to regulate the fraction to the heat exchanger is described by the 
following relations: 

For |t ! - T s | < 0. 1°F 

the fraction remains unchanged. Otherwise; 


f 2 - f x 

f - fi + 0.025 Tj - T s 

Ti - T 2 


and f is clamped between 0 and 0.9. 
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Where: f - new flow fraction bypassing heat exchanger 

f 1 « last flow fraction 
f 2 - flow fraction prior to f^ 

Tl - last cabin temperature, °F 
T 2 - cabin temperature prior to Tj, °F 
T s - set point cabin temperature, °F 

This is tantamount to a straight integral control technique with a 
varying integration gain constant. 

The effects of this control on cabin temperature and humidity are 
negated, however, by the logic In GPOLY2. GP0LY2 simply sets the 
cabin temperature to 70°F, the humidity to 42%, and the dew point 
to 46°F. This was done to give reasonable temperatures until a new 
control low can be developed. 

B.3.7.5 Water Tank Level Control 



Four water tanks are used for storage of clean and waste water; 
they are: 

(1) Urine and Wash Water Storage Tank 

(2) Clean Hygiene Water Storage Tank 

(3) Potable Water Storage Tank 

(4) Condensate Water Storage Tank 
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For each of these tanks, the calculations for the entering and 
exiting water flows are presented in the following paragraphs. 


B. 3. 7. 5.1 Urine and Wash Water Storage Tank 

The Urine and Wash Water Storage Tank receives water from crew 
urination, hand washing, and showering. The amount of water for 
each of these activities is considered to be dumped into the tank 
over a period of one time step which is currently 120 seconds or 
two minutes. Table B-5 gives the amount and time of day at which 
water is dumped into the tank for each of these activities. This 
table, of course, may be altered by the user by appropriate changes 
to the logic in GP0LY1. 

Flow will exit the tank only if the VCD unit is on. The VCD unit 
turns on whenever the tank level is greater than 30 percent full 
and turns off if the tank level falls below 27 percent full. When 
the VCD is on, it draws water at 2.3 lbm/hr for a three man unit 
and 2.73 lbm/hr for a six man unit. 

B.3.7.5.2 Clean Hygiene Water Storage Tank 

The Clean Hygiene Water Storage Tank receives water from the VCD 
unit at the rate processed by the VCD. Water is used from this 
tank for handwashing and showers and to supply the needs of the SPE 
and Multifiltration units. Accordingly, water for handwashing and 
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TABLE B-5 

SCHEDULE FOR DUMPING INTO URINE AND 
WASH WATER STORAGE TANK 


Water Dumped (Ibm) 


Time of 
Day 

8:06 AM 
9:00 
10:00 
11:00 
12:00 
1:00 PM 
2:00 
3:00 
4:00 
5:00 
6:00 
7:00 
8:00 
9:00 
10:00 
11:00 
12:00 
1:00 AM 
2:00 
3:00 
4:00 
5:00 
6:00 
7:00 
8:00 


3 Men 

Urine 

6 Men 

Hand- 

wash 

2.429 

3.754 

1.15 

2.429 

3.754 

1.15 

2.429 

3.754 

1.15 

1.15 

2.429 

3.754 

1.15 



1.15 

2.429 

3.754 

1.15 

1.15 

2.429 

3.754 

1.15 

2.429 

3.754 

1.15 

17.00 

26.28 

11.5 


Shower 


22.5 


22.5 
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showers is drawn from the tank according to the schedule in Table 
B-5. Again, all the water for a given time in the table is 
presumed to be drawn over one time step, i.e., two minutes. 


The water required by the SPE unit is given in Section B.3.1 by the 
relation for m^o.t* The water required for the multifiltration 
unit is 2.55 lbm/hr. However, water is drawn from the Hygiene Tank 

only if the condensate tank is unable to supply the multifiltration 
needs. 


B.3.7.5.3 Potable Water Storage Tank 


The potable water storage tank receives water from the multifiltra- 
tion unit. The multifiltration unit is supplied by a water pump 
which draws water from the condensate storage tank or the clean 
hygiene water storage tank if the condensate tank has insufficient 
water. The pump draws water at 2.55 lbm/hr. The multifiltration 
unit pump turns on when the tank level falls below 90% full and 
turns off if the level rises above 45%. 

Water is drawn from the tank for drinking and food preparation. 
Again, the amount of water drawn for each of these activities is 
considered to be drawn from the tank over a period of one time step 
which is currently two minutes. Table B-6 gives the amount and 
time of day for each of these activities. 
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TABLE B-6 

USAGE SCHEDULE FOR POTABLE WATER TANK 


Water Used (lbm) 


Time 
of Day 

Drinking 

3 Men 

6 Men 

Food Preparation 
3 Men 6 Men 

8:06 AM 
9:00 

2.794 

4.474 

1.185 2.370 

10:00 

11:00 

2.794 

4.474 


12:00 
1:00 PM 
2:00 

2.794 

4.474 

1.185 2.370 

3:00 




4:00 

5:00 

6:00 

2.794 

4.474 

1.185 2.370 

7:00 

8:00 

2.794 . 

4.474 


9:00 

10:00 

11:00 

12:00 

2.794 

4.474 

1.185 2.370 

1:00 AM 
2:00 




3:00 

4:00 

5:00 

6:00 

7:00 

8:00 

2.794 

4.474 


19.560 

31.320 

4.740 9.480 
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B.3.7.5.4 Condensate Water Storage Tank 


The condensate water storage tank receives condensate water from 
the condensing heat exchanger and the Sabatier reactor water 
separator. Water Is drawn from the condensate tank as required by 
the multif i ltration unit to supply the potable water needs. When 
water is drawn, it is drawn at 2.55 lbm/hr. 
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C.l Introduction 


This manual provides the model description document for a Space 
Station model which includes a habitat, laboratory, and four nodes. 
Only the air revitalization equipment is modelled; waste water 
management tanks and processing equipment are not included in the 
model. The principal pieces of equipment are: 


Function 


Subsystem Option Available 


CD 2 Removal 
032 Reduction 
O 2 Generation 
Trace Gas Removal 
Condensate Processing 


EDC, Molecular Sieve 
Bosch, Sabatier 
SPE, KCH 

Catalytic Oxidizer 

Plate-fin shuttle type heat exchanger 


Options are also available for hydrogen or 0C>2 bussing. 
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Figures C-l through C-14 present the schematics of the system as 
modeled using G189A. The following discusses how this system 
represents real hardware. Further descriptions of the system can 
be found in the User’s Manual [4]. 


In these schematics, extra lines, mixers, and splitters are 
inserted to provide the user with flexibility to select various 
options without having to rewrite the programs. These lines do not 
represent actual plumbing arrangements of a Space Station. For 
example, a duct does not exit the habitat then tee to two modes as 
shown in Figure C-l. In actuality, a node is adjacent to the 
habitat and a fan draws air directly from the habitat into the 
node. Therefore, these are functional schematics and do not 
represent actual plumbing. 
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FIGURE C-l 

SPACE STATION MXEL OVERVIEW 
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FIGURE C-2 
SPACE STATION NODES 
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FIGURE C-3 

OVERVIBV OF HABITAT ARS 
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FIGURE C-4 

OVERVIEW OF IABORATORY ARS 
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FIGURE C-13 

HABITAT CATALYTIC OXYDIZERS 
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FIGURE C-14 

LABORATORY CATALYTIC OXIDIZERS 
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Modellin g of Components 

The following sections discuss the analytical 
components available in this Space Station model, 
have been described elsewhere. The following lists 
and where they have been discussed. 


Subsystem 


Reference 


SPE 

Catalytic Oxidizer 

Sabatier 

EDC 


Appendix B 
Appendix B 
Appendix B 
Appendix B 


The following subsystems are discussed below: 


Molecular Sieve 

Bosch 

KCH 

Plate Fin Heat Exchanger 
Control Logic 


models of the 
S ome of these 
the subsystems 


C— 17 


SVHSER 10638 



UNITED 

TECHNOLOGIES 

GmME3T@Kl 


C.3.1 


Molecular Sieve 

Molecular Sieve removes carbon dioxide from incoming air. S ome 
constants are: 


*min 

t nms 

t msg 

cfm 

A Pfan 


^ Hsg 

L %s 

CH^ 


- half cycle time — 60 minutes 

= min. pressure of desorbing mole sieve bed = 1.0 psia 

- max. temperature of desorbing mole sieve bed = 360°F 

- max. temperature of desorbing silica gel bed - 180°F 
-fan cfm = 23.25 

- 11 in H 2 0 

-fan efficiency =0.35 

- OC >2 removal efficiency - 0.65 

= enthalpy change for adsorbing silica gel bed = 1400 
Btu/lb H 2 0 

- enthalpy change for adsorbing mol sieve bed - 1800 
Btu/lb GC^ 

- compressor cfm - 1.05 


In order to understand the analytical description, the schematic of 
the molecular sieve subsystem shown in Figure C-15 should be 
reviewed . 
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First, the temperture leaving the fan is computed by the following: 


where: T 2 

Tl 
Qfan 


t 2 


Qfan 
Tl 


m a c p 

temperatures leaving fan, F 


temperature entering fan, F 
fan power, Btu/lb 


cfm 4 Pfan 

QFAN 3.413 

8.5 1 


m a - air inletmass flow to fan, lbm/hr 
Cp ■ specific heat of air, Btu/lbm-F 


Next, the temperature of the process air leaving the adsorbing 
silica gel bed is calculated from the following: 


LO Ahjg 

t 3MAX - t 2 + 

C P 

t 3MAX ~ t 3 

m sg ■ 

600 - t h 

T 3 ■ T 3 + m sg A t 

where: 00 - absolute humidity of air entering molecular sieve 

T 3 max * maximum temperature of air leaving the adsorbing 
silica gel bed, F 

m sg = change of air temperature leving silica gel bed 
with time, F/sec 
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At = simulation time step, sec 
th - time into half cycle, sec 


After 10 minutes into the half cycle, the temperature of the air 
leaving the silica gel bed is limited to T 3 max . 


The water removed by the silica gel bed is: 


m H20 ■ m H 20 + ( M H20,in) ^ * 

where: Mh 20 ■ mass of water removed by silica gel bed during 

half cycle, lbm 

Mh 20 in “ flow °f va P or an d entrained liquid entering 
system, Ibm/hr 

Thus the assumption is that all the moisture entering the system is 
adsorbed onto the silica gel bed. Now, the temperature and mass 
flow of cool dry air leaving the heat exchanger is: 

T4 - T c00 i + 5°F 

where: T coo i - temperature of coolant water entering HX, F 

T 4 « temperature of air leaving HX, F 
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The heat transferred by the heat exchanger to the water is: 


Qhx = m d c p ( t 3 ” t 4) 

where: - mass flow of dry air which entered molecular sieve 

subsystem, Ibm/hr 

The air flows next to the CO2 adsorbing molecular sieve bed. The 
temperature of the air leaving that bed is given by the following 
calculations: 


m C02,a a m C02,i £ C02 

M C02,a ■ M C02,a + m C02,A A t 

p C02,e - p C02,i (1- £c02> 

m C02,e " m C02 , 1 (1 ~ £ C02) 

M C02,a Ah ms 
t 5,MIN - t 4 + 


m d c p 


T 5,min " t 5 

%s * 

1800 - t h 


t 5 ■ t 5 + %s At 
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where: mco 2 ,i - inlet c02 mass fl0W int0 molecular sieve bed ’ pph 

£ C02 “ CO 2 removal efficiency of molecular sieve bed 
m C02,a ■ rate of co 2 adsorption, lbm/hr 
M C02,a * total mass of co 2 present on bed, lbm 
Pcq 2 i * partial pressure of CO 2 in inlet air, psia 
p C02,e ■ partial pressure of CO 2 in molecular sieve exit 
bed, psia 

m C02,e * exit co 2 mass flow out of molecular sieve bed, pph 
T 5 ,min - minimum temp, exiting molecular sieve bed, F 

After 1,800 seconds into the half cycle, the exit temperature 
cannot be less than T 5 >nV j n . 

The process air then flows to the desorbing silica gel bed. The 
temperature of the bed rises for the first 17 minutes, peaks, then 
falls for 17 minutes until the temperature reaches a minimum 
T 5 ,min* The equations are: 

For the first 17 minutes: 


M sgd 


T MSG " t 6 
1020 - t h 


T 6 - T 6 + M sgd At 

where: ts - temperature leaving, desorbing silica gel bed, F 

M S gd ■ temperature change with time of airleaving desorbing 
silica gel bed, F/sec 
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This temperature Tg cannot exceed T mS g. 

For the next 17 minutes: 

T 5 ,min - Tg 

M sgd * 

2040 - t h 

^6 " ^6 + Mgg^j At 

For these 17 minutes, the air temperature is limited to a minimum of 
T 5,min* From 34 minutes into the half cycle to the end of the half 
cycle, the temperature exiting the desorbing silica gel bed is set to 
T 5,min* 

The properties of the air leaving this silica gel bed are: 

P H20 - p SAT( t 6) 


mw H20 P H20 

m H20 m d 

MW d P T " P H20 

M sgd ■ M sgd + m HS0 At 

where: P HS o - partial pressure of water vapor in air, psia 

m H20 ■ mass flow of water leaving silica gel bed, pph 
M S gd - mass of water desorbed from bed this half cycle, lbm 
Of course, the mass desorbed is limited to the mass adsorbed from the 
previous cycle. 
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The electrical power to the desorbing silica gel bed is given by: 


For tft < 240 seconds 
th > 240 seconds 


W htr - 0 watts 
W htr « 657 watts 


The average heat given up by the desorbing silica gel bed to the 
cabin is: 


46 

Qsg ” — ( T msg “ T cab) 
y 290 

where: T^j^ * cabin air temperature, F 

The final calculation is that for the mass of C0 2 desorbed from the 
molecular sieve bed. The following equations are used for t h > 480 

seconds. 


R (T 4 + 460) 

V J 

144 P min MW C 02 

( P8 \ 0 * 769 cfm c 

M C02,d = C 1 * 01 + °* 01 (" - 1 Y7~ 

V^miny ' 

where: R - Universal gas constant = 1545 ft-lbf/mole R 

MWco 2 * Molecular weight of C0 2 = 44 lbm/mole 
P m ^ n > Minimum pressure of desorbing molecular sieve, psia 
Pg - Accumulator pressure, psia 
cfm c - Compressor cfm c 
m C02,d “ Mass flow of c0 2 desorbed, lbm/hr 
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For th < 480 seconds, mco2,d * 0. 

The total carbon dioxide desorbed during the present half cycle is 
given by: 


M C02,d “ M C0 2,d + mC02,d 4t 

Of course, the total amount desorbed is limited to the amount of 
C0£ that was originally present on the bed. 

C.3.2 Bosch 


The Bosch is a process that reduces carbon dioxide and hydrogen to 
carbon and water while giving off heat. The model used here is 
based on the one described in the G189A Manual [2]. A functional 
schematic of the process is shown in Figure C-16. The chemical 
reaction is described by: 

2H 2 + C0 2 > C + 2H 2 0 

A listing of the program is provided in the User's Manual 
Appendices [4]. 

First the molecular weight of bone dry condenser exit gas and the 
inlet molar flows of H 2 and C0 2 are computed. 
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FIGURE C— 16 

BOSCH PROCESS SUBSYSTEM SCHEMATIC 
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n H2, i - m H 2 , 1 /2.016 
n C02,i * m C02 , 1 ^ 44 .011 

MW BDG - 16.043 Y CH4 + 44.011 Y C0 2 + 2.016 Y H2 + 28.011 Y co 

where: m H2,i “ Mass flow of H 2 entering Bosch system, pph 

n H2,1 * molar flow of H 2 entering Bosch system, moles/hr 

m C02 , i * mass flow of C0 2 entering Bosch system, pph 

n C02, i a molar flow of C0 2 entering Bosch system, pph 
y CH 4 ■ mole fraction of methane in dry cond. exit gas = 
.235 

Y G o 2 = mole fraction of C0 2 in dry cond. exit gas = 

.163 

Y h2 - mole fraction of H 2 in dry cond. exit gas = .327 

Yqq - mole fraction of CO In dry cond. exit gas - .275 

MW B qq * molecular weight of bone dry cond. exit gas 

The amount of carbon processed and water produced for a 
quasi-equilibrium assumption Is given by the following which is 
based on the reaction formula: 

if C0 2 1 imi ting: 

m C,p “ m C02 , i * MW C 
m H20,p ■ 2 m C02 , i * mw H 20 + m H20,i 
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If H2 1 imiting: 


m c,p = m H2,i MW c /2 
m H20,p * m H2, i mw H 20 + m H20,i 

where: MW C - molecular weight of carbon * 12.011 Ibm/mole 

MWh 20 ■ molecular weight of water * 18.016 lbm/mole 
M H20,i M mass flow va P° r and entrained liquid entering 
Bosch system, lbm/hr 

m c> p » mass rte of carbon production, lbm/hr 
m H20,p ■ mass rate water production, lbm/hr 

Therefore, flow out of condenser Is: 

m 3 « m r - m CjP - m H 2 o,p 

where: m r * recycle flow rate on more specifically the mass 

flow at the compressor « 6.80 lbm/hr 

The temperature at the condenser exit is iterated upon. Its 
initial value is assumed to be 20°F hotter than the inlet coolant 
temperture: 


t 3 * T cool + 20 
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The iteration begins with the calculationof the flow rate of bine 
dry condenser exit gas 


m BDG 


P H20 - p sat ( T 3) 


1 + 


m 3 

mw H20 p H20 

mw BDG p C0ND " P H20 

where: Ph 20 * partial pressure of vapor leaving condenser, psia 

Pcond * total pressure of recycle gas in condenser = 

16.9 psia 

mBDG * mass flow bone dr y 9 as leaving condenser, pph 


Accordingly, the flows out of the condenser and compressor are: 


m BDG 
n BDG “ 

mwbdg 

P H20 

m H20, 3 " n BDG mw H20 

p cond 

m H20, 4 - m H20, 3 + m H20, 2 
m H2, 3 - n BDG y H2 mw H 2 
m H2, 4 - m H2, 3 + m H2, 4 
m C02, 3 * m BDG y C 02 
m C02 , 4 - m C02 , 3 + m C02, 4 
m CH4 , 3 - n BDG y CH4 mw CH 4 
m CH4, 4 “ m CH4 , 3 

m C0, 3 M n BDG Y C0 MN C0 

m C0, 4 * mCO, 3 
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The properties of the flow exiting the compressor are: 


Cp 4 = (3.42 mn 2,4 + 0.21 mQQ 2,4 + 0.55 mgn4,4 + 0.25 nigg^ 

+ 0.49 m H 20,4 + 0.22 mg 2 + 0.299 m N 2 > / m r 

MW 4 =■ m r / (mH2,4/2.016 + m C 02 , 4/44. Oil + ""CM, 4/16. 043 
+ m C0, 4/28. Oil + m H20,4/18.06 + m 02/32 
+ m N2/28.088) 


y~A 


where: 


1,987 



Cp 4 - Specific heat of gas at compressor exit, BTU/lbm-F 
MW 4 - Molecular weight of gat at compressive exit, 
lbm/mole 

7*4 * Ratio of constant pressure to constant volume 
specific heats 

mo 2 ■ Mass flow of oxygen In recycle gas - 0.0 lbm/hr 
m ^2 - Mass flow of nitrogen in recycle gas « 0.0 lbm/hr 
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The mixed temperature of inlet gases and the recycle gases exiting 
condenser is given by: 


T 6 = C( m BDG 1 Cp BDG + °*49 m H2 o,3) T 3 + °* 21 m C02,i + 3 * 42 m H2,i 
+ 0.49 m H2 o,iv + 1-0 m H20,ie] / ( m r C P4* 

Cp BDG - C 3 * 42 m H2,3 + °- 21 m C02 , 3 + °- 55 m CH4,3 + 0.25 m C0,3 
+ 0.22 m 02 , + 0.299 n)N 2 )/m BD 6 


where: 


Tg - Mixed temperature of gases entering, compressor, F 
m H20, iv ■ Mass flow of vapor entering system, pph 
m H20, ie ■ Mass flow of entrained liquid entering system, pph 
Cp BDG * Specific heat of dry gas exiting condenser, Btu/lbm-F 


The energy required to raise the pressure of the recycle flow gases 
across the compressor is given by: 



r n - 1 

. 

Qc * m r R (T G + 460) r 4 

(-=--) V. 

-1 


_' p cond' 



778 MW 4 (y 4 -l) fy'a 

where: R - Universal gas constant - 1545 ft-lbf/mole-R 

P r - Reactor pressure = 24.3 psia 
a ■ Aerodynamic efficiency of compressor - 1.0 
Q c = Compressor energy into recycle gas, Btu/hr 
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The total power consumed by the compressor is given by: 


Qc 

W c - r 

e c 3.41 


where: W c - Power to run compressor, watts 

£ c - Motor efficiency of compressor 


Therefore, the heat lost to the ambient is: 


Q a - 3.41 W c - Q c 


The heat of reaction in the reactor is: 

Q r - 973 x mco2,i 

Essentially, it is assumed that the input gases are stoichiometric 
in ratio which then causes the mole fractions of the gases to 
remain constant in the recycle loop. The temperature out of the 
compressor is: 


T 4 


T6 + 


Qc 

m r Cp 4 
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The temperature out of the reactor is: 


m r Cp4 

where: Q|_r = Heat loss from reactor - 0.0 Btu/hr 

The temperature into the condenser is: 

t 2 - T 1 " ^HX ( T 1 - t 4) 

where: ■ Heat exchanger efficiency - 0.85 

Finally, the temperature of the gases leaving the condenser is: 

t 3 - t 2 ~ ^c ( t 2 ~ T cool ) 

where: € c - Condenser efficiency - 0.90 

This T 3 is compared with the original gressed T 3 . If they agree 
within 0.3 F, the iteration Is complete; otherwise, this T 3 is 
tried as the next gress. 

C.3.3 Static Feed Water Vapor Electrolysis (KOH) 

The Static Feed Water Vapor Electrolysis subsystem uses KOH as the 
medium for electrolysis to produce oxygen. This model is based on 
the ELCELL subroutine in G189A [2]. The following assumptions are 
used: 

(1) The unit is isothermal within the cells. 
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(2) Product gas streams exit at the prevailing cell temperature. 

(3) The cell Faradaic efficiency is 100%; the cells are voltage 
inefficient only. 

(4) All cells are connected -in series, i.e., the same input 
current passes through each cell. 

(5) The thermoneutral voltage is 1.48 volts; the voltage 
efficiency is equal to the thermoneutral voltage divided by 
the actual voltage. 

The gas production is calculated as follows: 

J - I/A 

m H20,e * N c 1/1350 
m H2,p “ M H20,e/ MW H20 
n 02,p * n H2,p /2 
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where: I - Input current to cells, amps 

A - Area of one cell, ft2 
N c - Number of cells 

m H20, e * Mass flow 02 consumed by electrolysis lbm/hr 
J = Cell current density, amps/ft2 
n H2, p ■ Molar flow of hydrogen produced, moles/hr 
n 02, p “ Molar flow of oxygen produced, moles/hr 


The details of the analytical method are described fully in the 
G189A manual for the subroutine ELCELL. However, the calculation 
of the energy required for electrolysis is computed according to 
the following equations and tables: 


First, the cell voltage at 150 amps/sq. ft. current density for any 
cell temperature is determined by interpolation of the following 
table: 


'C 

ill 

»0 

(Volts) 

110 

1.660 

120 

1.620 

130 

1.575 

140 

1.550 

150 

1.530 

160 

1.500 

170 

1.475 

180 

1.465 

190 

1.450 


C-36 



SVHSER 10638 



UNITED 

TECHNOLOGIES 

GMKfflOTOK] 


For cell temperatures less than 110 F, the voltage is set to 1.70 
volts. 


Next, the effect of a current density different from 150 amps/sq. 
ft. is found by interpolating the following table: 


J Delta V/ Delta J 

(amps/ft) (Volts/ASF) 

0.01460 
0.00060 
0.00055 
0.00050 
0.00050 
0.00045 
0.00045 

Then the following is used to calculate the cell voltage: 


0 

100 

200 

300 

400 

500 

600 


V c - V 0 + _V (J-150) 


The cell efficiency is: 


Y) = IAQ ioo 
V c 


and is limited to a peak value of 99%. 


Lastly, the energy required for electrolysis is given by: 


Q - N c I V c 3.413 


C.3.4 Plate Fin Condensing Heat Exchanger 


This subroutine models the performance for a plate fin condensing 
heat exchanger. Basically, the program Iterates on the condenser 
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C.3.4 Plate Fin Condensing Heat Exchanger (Continued) 

outlet temperature until the outlet temperature does not change 
from iteration to iteration. 

First, the inlet dew point is calculated: 

m v i 

u> i „ --1- 
m o,i 

w i p t,i 

P v i 1 

s W i + 0.622 
T dp,i * T sat ( p v,i) 

where: m v j * Mass flow of inlet vapor, pph 

m da,i * Mass flow of dry air, pph 
^ 1 » Inlet absolute humidity 
Pjj - Inlet air total pressure, psia 
p vi ■ Inlet air partial pressure of water, psia 

The initial guess at the exit temperature is one quarter of the way 
from the coolant inlet temperature to the air inlet temperature: 

T a,e ■ T c,i + ( T c,i + T a,i) 

From the exit temperature, the outlet absolute humidity and vapor 
pressure are calculated: 
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p v,e = p sat ( T a,e) 

<^ P - 0.622 p v,e 

P T - p v,e 

The heat loads are then: 

Ql * Mda,i ( w i - w e) h fg 
Qs - m a , 1 c p,a ( T a,1 ~ T a,e) 

Qt ■ Qs + Ql 

where: hfg - Heat of vaporization, Btu/lbm 

Ql_ - Latent heat, Btu/hr 
Q s - Sensible heat, Btu/hr 
Qf - Total heat loss, Btu/hr 
m a j * Total mass flow of air Into Hx, lbm/hr 

The coolant exit temperature Is: 

Qt 


where: m c - Coolant flow, pph 

Cp.c - Specific heat of water, Btu/lbm-F 

From here, various properties of the air and coolant are calculated 
as a function of the average temperature. The properties and 
nondimensional numbers calculated are: 
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m - Viscocity 

K - Thermal conductivity 

P r - Prandtl number 

R e - Reynolds number 

C 0 i - Colburn factor 

C m od ■ Modified Colburn factor 


The conditions are shown in the listing presented in the User's 
Manual [4]. Also shown there are various geometric dimensions such 
as fin height. 


The film coefficient is chosen as the maximum of the following two 
equations: 


h - 3.65 K/D h 

h - C m0 d C p G P r --6667 

Overall fin efficiency is then computed from the following fin 
equations: 



c - 2- 
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The effective UA becomes: 

UA - ^ hA 

The total effective dry UA is calculated as 

1 

UA dry “ ■' j' 

UAc + UAj 
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UA C - Cold side effective UA, Btu/hr-F 
UA a - Hot or air side effective UA, Btu/hr-F 


The pinch point temperatures for the air side and coolant side are 
given by: 


PP.a 


/UAa s 

* T «p.1 - T c,e] C p .c i"c + m a t i Cp, a T a>1 


m, 


'a,i c p,a + m c C p#( 


UA, 


UA r 


The pinch point temperatures for the air and the coolant are those 

which occur at the location where the wall temperature equals the 
inlet dew point temperature: 


T PP.c " T dp,i ~ UA a (T pp>a - T^pj) 

UA C 

These above equations can be derived from the simultaneous solution 
of the following two energy balances: 


UA a (Tpp.a - T dPi1 ) . UA( . (T<)pj1 _ Tpp c) 
m a,1 Up. a (Tpp.a - T a>i ) , m c Cp>c (Tpp>c _ Tce) 

The wet side and dry side log mean temperature differences are: 


C-42 




UA d,tot * UA dry 
C.3.5 Control 

The control of the Space Station model is done in subroutine 
GP0LY1. The following paragraphs describe the control laws used to 
control : 
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(a) Oxygen partial pressure 

(b) Total cabin pressure 

(c) Oxygen accumulator pressure 

(d) Temperature control 

(e) C02 accumulator exit flow control 

(f) H2-C02 mix to C02 reduction unit 

Oxygen partial pressure is controlled by using the same technique 
as described in the original ESCM Model Description Document [1]. 

The controller maintains oxygen partial pressure between 3.09 and 
3.23 psia. 

Total cabin pressure Is maintained by the addition of nitrogen. 
The controller admits nitrogen to bring the pressure up to 14.813 
psia only after the oxygen pressure is above 3.09 psia. The 

original Model Description Document [1] should be seen for more 
detai 1 . 

The oxygen accumulator pressure is maintained by addition of oxygen 
from the oxygen generators as oxygen is drawn from the accumulators 
to maintain cabin 0 2 partial pressure. The current to the oxygen 

generators is regulated to each of the generators according to the 
following law: 
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For Pq 2 2. 1050 
For Pq 2 i 950 
For 950 P 0 2 < 1050 

As the current is increased, 
electrolysis unit increases. 
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1 * 0.9 I nom 
I - 1.1 I n0 m 
1 * Inom 

the oxygen generation by the 


For temperature control, the bypass flow around the condensing heat 
exchanger is regulated. A proportional plus Integral scheme is 
used: 


R - 0.05 E + JO. 0001 E dt 

where: R - Fraction of flow to bypass Hx. 

E - Temperature error - T se t - T ac t 

As more flow bypasses the Hx, the mixed flow after the Hx is 
hotter. 

Carbon dioxide removed from the air is sent to an accumulator. 
Each C02 removal unit has its own accumulator in a nonbussed 
system. The flow out of the accumulator is regulated to maintain 
pressure in the accumulator above 21 psia. Essentially, the 
average C02 removal rate during the past molecular sieve cycle is 
used as the flow out of the accumulator for the present cycle. 
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Lastly, the H2-C02 mixture into the C02 Reduction unit is regulated 
to be stoichiometric by venting excess hydrogen after being 
produced by the oxygen generators. The required H2 flow and the 
amount to be vented is given by: 


MW H2 

m H2 “ 2 m C02 

MWqq2 

m 1 - m H2 
R - m^ - 

«»1 

where: mco2 “ Flow of C02 to reduction, pph 

m^2 ■ Stoichiometric flow of H2 to reduction, pph 
m^ * Inlet flow to splitter, pph 
R « Fraction of H2 inlet flow to be vented 
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